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Abstract 

In this thesis exact results on 0(a 2 ) single bremsstrahlung corrections to low 
angle Bhabha scattering at LEP/SLC energies are given. The calculation rep- 
resents the last outstanding theoretical second order subleading electroweak 
contribution for that process, needed to determine the experimental luminos- 
ity at second generation LEP detectors below the 0.1% precision threshold. 
The exact, fully differential result is obtained by employing analytical as well 
as computer-algebraic methods and includes terms up to 0(0.05%) relative to 
the Born cross section. The initial output of over 20, 000 terms could be re- 
duced to 90, only 18 of which are shown to be numerically relevant and for 
which a simple logarithmic ansatz is derived, that is in remarkable agreement 
with the complete answer. Strong consistency checks are performed, including 
Ward-Takahashi identities and tests on the right infrared limit according to the 
Yennie, Frautschi and Suura program. Monte Carlo results for the integrated 
cross section are compared with existing calculations in the leading logarithmic 
approximation for a chosen set of experimental cuts. The size of the missing 
subleading terms is found to be small but non negligible in the context of setting 
stringent limits on Standard Model predictions and thus its realm of validity. 
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Chapter 1 

Introduction 



High precision measurements at the LEP/SLC colliders have made tremendous 
improvements over the last few years. At LEP, for instance, first generation 
detectors have measured the absolute luminosity with an accuracy of 0.3 — 0.5% 
||. Present detectors have reached a precision below the 1pm threshold ||. 
This level of precision opens up a wide range of new stringent constraints that 
can be specified not only for physics beyond the Standard Model, but also for 
the Higgs sector and other parameters of the electroweak theory (T^]. While 
LEP II is now on line to probe the W + , W~ aspect of the GSW- model, there 
is still a lot of interest in the LEP / energy regime of the intermediary neutral 
vector boson. 

In general, the Z line shape is described by three parameters: Mz, I 'z and cr° 
||. The luminosity normalizes the line shape cross section and it thus directly 
influences the accuracy of the <t° measurement. In order to determine the lu- 
minosity L, one needs to count the number of events N for a process in which 
the cross section a is "known": 

L = N - (1.1) 
a 

The "known" process at LEP/SLC is low angle Bhabha scattering p8| . More 
than 99% of this cross section is due to t-channel photon exchange J44| and 
can in principle be calculated in QED with arbitrary precision. The remaining 
part of the cross section is due to s-channel Z, t-channcl 7 interference, and its 
relative contribution decreases with decreasing scattering angle H . 
The overall accuracy of the "known" cross section is very important for high 
precision measurements of Standard Model parameters since a relative luminos- 
ity error ^7 = 10 -3 , for example, changes oq by 42p6 and N u , the number of 
massless neutrino generations, by 0.0075 ||. Also T e = T(Z — > e + + e~) is 
directly affected by the precision of L since |3^] 

r.' - ^ , (1.2) 

1Z7T 
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CHAPTER 1. INTRODUCTION 



which is one of the main sources of our knowledge of the electroweak mixing 
angle |52| . It is is thus necessary to have theoretical calculations that can match 
the experimental accuracy below the lpm mark. As can be seen from Refs. 
|35| , |39| , pic| ] , this has not yet been achieved, although initial state s-channel 0(a 2 ) 
corrections, including two loop effects have been calculated The result of 
Ref. Q for one virtual and one hard photon, however, is not suited for a Monte 
Carlo implementation since it was calculated in an "inclusive" way, i.e. it is 
not differential in the photon angles, which prohibits any adaptation to a given 
detector geometry. This is reflected in the form that the aspect of the theoretical 
uncertainty that contributes most to the present 1.6pm precision level is given 
by the missing parts of the second order subleading terms originating from 
the one photon bremsstrahlung process [ p5| . In this thesis, the exact O(a) 
virtual corrections to the single hard bremsstrahlung cross section at = Mz 
are calculated up to the precision required in the above context, i.e. up to 
0(< 0.05%). This means, for instance, that two photon exchange contributions 
can be neglected |34| as well as t-channel Z exchange diagrams [[35| . The result 
is calculated in differential form and Monte Carlo results then give the 
integrated cross section that is relevant for the theoretical prediction of the low 
angle Bhabha scattering cross section, needed to determine the luminosity L. 
This work is structured as follows: 

Chapter § describes the conventions and defines the notation used in this the- 
sis. A brief section is devoted to scalar integrals and how to approach complex 
loop integrals in the context of using an algebraic manipulation language. It 
also gives an expression for the single bremsstrahlung tree level amplitude A T ^- , 
which will prove to be pivotal for the final form of the exact differential result 
given in chapter ^. 

In chapter ||]the details of the calculation of the 0(a) virtual corrections to the 
single hard bremsstrahlung cross section are presented. The nature of the prob- 
lem of evaluating higher order Feynman amplitudes brings with it a somewhat 
unavoidable technical aspect in the various sections. This chapter contains the 
new results presented in this work in conjunction with the appendices. 
A summary of the various expressions of chapter ^ is given in chapter ^ as the 
complete differential result. First numerical results are presented and discussed. 
Chapters |B| and |6| are devoted to demonstrating that the presented exact 
differential result satisfies various internal consistency checks and has all the 
properties expected from general theoretical arguments. 

The total cross section resulting from the differential expression in chapter ^ 
is obtained by means of a Monte Carlo phase space integration and numerical 
comparisons with other published leading log (LL) type calculations are pre- 
sented in chapter for cuts similar to those used at the luminosity detector 
SICAL at ALEPH in LEP [j36f . Also, a new approximation is derived that shows 
excellent agreement with the exact result and the overall size of the subleading 
terms is discussed. 

Finally, chapter ||| contains a summary of the results obtained in this thesis in 
the context of the present status of high precision radiative corrections and also 
concluding remarks. 



Chapter 2 

Basic Definitions and 
Identities 



In this chapter the basic notations and definitions will be discussed as well as 
some identities that follow from the high energy limit. If nothing else is stated 
explicitly, the conventions of Bjorken and Drell are used throughout this 
work. Since the overall objective of this thesis is a high precision calculation of 
0(a 2 ) radiative corrections to a high energy process, it will prove to be very 
convenient to make use of the so called high energy limit, where terms of order 

2 

can be neglected. As a consequence of rendering the electrons (positrons) 
massless, their helicity will be conserved. This then leads to the very useful 
concept of introducing fermion helicity states |^| which are defined here by the 
following Dirac notation 



«a(p) 



v-x(P) 
v-x(P) 



(2.1) 
(2.2) 



Eqs. 2.1, 2.1 hold for arbitrary four-momentum P with P 2 = and u, v defined 
as m Ref. [§. The normalization is fixed by demanding that 



(P,A| 7/1 |P,A) 



2P„ 



(2.3) 



For two arbitrary massless four- momenta, P and Q, the following equations 
hold: 



{P,X\Q,X) = 0; (P,-A|P,A) - 0;P|P,A> - (2.4) 

(P,-A|Q,A) = -<Q,-A|P,A);P = |P, A)<P, A| + \P,-X)(P, -A| (2.5) 

(P,-X\Q,X)(Q,X\P,-X) = 2PQ = (P + Q) 2 (2.6) 

Furthermore, for arbitrary four-momenta kf. 
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(P,\\fi 1 ..4 n \Q,X) = (Q,-A|^„...^i|P,-A) ,n odd (2.7) 
(P,-A|#i...#„|Q,A) = -(Q,-A|^ n ...^i|P,A) ,neven (2.8) 



An important identity that will be used extensively is the Fierz-identity 19 
and holds for arbitrary spinors \ A,X) , \B,X') , \C,p) and \D, //) : 



(A A| 7 , |B, A) (C, M | 7 V | A M) = 2 (A, A|X, -A) (Y, -A|B, A) , (2.9) 

where 

It is sometimes convenient to simply write 

(P,Q) X = (P,-X\Q,X) (2.11) 
The "magic" polarization vector of Xu, Zhang and Chang is given by : 

e^K P ) = — _ , (2.12) 

where K is the photon four-momentum, p its polarization and h a reference 
momentum with h 2 — 0. Choosing h amounts to a gauge choice and this 
freedom can be exploited to greatly simplify calculations with bremsstrahlung 
effects |ll]]. If the incoming electron has a four-momentum P, the outgoing one 
a four- momentum Q and if its helicity is denoted by A, the "magic" choice in 
the present calculation turns out to be : 

For completeness, there are three other helicity dependent four-momenta used 
in this work, h', h and h', where analogously for the positron line 

h ' = {q> %Z- p ^ 

h and h' denote the helicity flipped choices of h and h! respectively. 

It is also helpful to give ^ explicitly, since it is the object appearing in expressions 

for Feynman amplitudes : 

(\K,-p)(h,-p\ + \h,p)(K,p\) (2.15) 



(h,K) p 

As an application, the "magic" properties of are used in proving the following 
important identity for A = p: 
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= PQ(Q',\'\t\P',\')(Q,\\fi\P,\) 

Qe {Q',X'\P\P',X'){Q,X\$\P,X) + 
PKQe (Q',X'\^\P',X') (Q,X\^\P,X) 



(2.16) 



Eq. 2.16 can be proved by observing that 



PK(Q,X\^\P,X) = P»{Q,X\fi\P,X) - ~{Q,X\P^$\P,X) (2.17) 



PQ(Q,X\fi\P,X) = ~{Q,X\PQfi\P,X) 
Thus, it remains to be shown that 



(2.18) 



= (Q',X'\t\P',X')(Q,X\PQfi\P,X) - 

Qe(Q',X'\ 7fl \P',X') (Q,X\P^fi\P,X) , 



(2.19) 



which follows from the transversality condition eK = and the "magic" prop- 
erties 



f\P,\) = & {$4} = ,X = p 
An analogous identity holds for A = — p and reads: 



(2.20) 







PQ{Q',Xy\P',X'){Q,X\fi\P,X) - 
Pe (Q',X'\(#\P',X')(Q,X\fi\P,X) + 
QK Pe (Q', A'| 7m \P', A') (Q, A| 7" \P, A) 



(2.21) 



Eq. 2.21 can be directly obtained from Eq. 2. 16| through the crossing rule 



P -Q , A 



-A 



(2.22) 



The usefulness of the helicity dependent spinor notations 2.13, 2.14 can be 



exemplified by stating the identities |2.16| , [2.2l| through a helicity independent 
notation: 



= hh{Q',X'\<t\P',X'){Q,X\fi\P,X) - 
he (Q , ,A / |^|P' ! A')(Q,A|^|P,A) + 
hKhe (Q',A'| 7ai |P',A') (Q,X\^\P,X) 



(2.23) 



Note in particular that the polarization vector e M in 2.12 is invariant under the 
transformation 2.22, which relates initial to final state electron-line radiation 
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and thus insures that the same polarization vector be used for one line in each 
helicity case. 

The following definitions will also b e ext ensi vely u sed throughout this work, 
where the four- moment a are those of |2.13 and 2.14: 



a = 


2PK 


= -{P-Kf 


(2.24) 


(3 = 


2QK 


= (Q + Kf 


(2.25) 


s = 


2PP' 


= {P + P'f 


(2.26) 


s' = 


2QQ' 


= (Q + Q'f 


(2.27) 


t = 


-2PQ 


= {P-Qf 


(2.28) 


t' = 


-2P'Q' 


= (P'-Q') 2 


(2.29) 



2 i 

It should be mentioned that the fine structure constant a = |- = will not 
appear explicitly in equations in this work unless specifically mentioned in order 
to avoid confusing notation. 



The t-channel process, as shown in Fig. 2.1, also allows for Z-boson exchange in 



the present calculation, mainly to be able to incorporate s-channel Z, t-channel 
7 interference and, at a later stage, to possibly allow for wide-angle applications. 
The propagator G\ t \(t') is then given by ]2S|| : 



q (f) 1 [(l-A)-4 S m 2 fl4(l-A')-4 S m\] 
A ' A ' 1 ' t' 4sin 2 20 w (t' - M z 2 ) ' 

where 9 W denotes the weak mixing angle and Mz the mass of the intermediary 
Z-boson |^0| . In the s-channel a correction coming from the imaginary part of 
the self-energy correction |||] needs to be taken into account and leads to the 
replacement rule Mz 2 — ► Mz 2 — * s if~j where Tz is the width of the Z-boson. 



2.1 Scalar Integrals 

All results that are given below involving loop integrals are expressed in terms 
of so called scalar integrals, i.e. four-dimensional integrals over Minkowski space 
with no tensor-structure in the numerators. The analytic solutions of integrals 
of this type have long been known in the literature |l7| and recently, a new 
way of numerically implementing these results in a more stable way has been 
proposed Jl3| . The results of Ref. are those that have been used in this 
work to calculate scalar integrals by means of a numerical package called //, 
as well as their reduction scheme, which decomposes tensor integrals into a 
linear combination of the easier scalar integrals. Wherever possible, the package 
results have been compared to known analytic results in certain limits JIt], [5l) 
and have shown very good agreement. The reduction scheme, implemented 
in the algebraic computer-language FORM has also passed all tests that 
were applied to it with relatively simple test-cases. It is essential for a numerical 



2.2. COMPLEX LOOP INTEGRALS AND BILINEAR COVARIANTS 13 



analysis that the algorithms, especially for the more complicated scalar integrals 
with three and four denominators, are reliable and indicate situations where 
the required cancelations become too involved to be handled within a given 
level of accuracy. This is one of the most important advantages the //-package 
possesses, where roughly 80% of the programmed code (« 50000 lines) is devoted 
to just that flllfl . The scalar integrals used in the present calculation are defined 
as follows: 



rk _ i r ifd n i 



iir 2 J (I 2 — mi 2 )((I + pi) 2 — m^)...((l + pi + ... + Pk-i) 2 - mk 2 ) + ie 

(2.31) 

Throughout this work only scalar integrals up to k = 4 will be necessary and 
are going to be denoted by 



^o 1 = 
Xi = 

Xq = Co{pi,p2 2 ,(pi+P2) 2 ,m 1 ,m 2 ,m 3 



Ao(mi) 

B (pi,mi,m 2 ) 



D (pi- 



2 2 2 
P2 ,P3 ,Pi ; 



(2.32) 
(2.33) 
(2.34) 

(P3+P4) 2 ,(pi + P4) 2 ,mi,m 2l m 3 ,m4) (2.35) 



It should be noted that, as suggested by the notation in Eq. 2.31, the scalar 



integrals 2.32, 2.33] are defined through the n-dimensional regularization method 



27, 21 . The choice of the scale /i should not affect any result. In order to write 
Eq. 2.31 in an easy to recognize manner for arbitrary k, all momenta pi are 
assumed to be incoming in a Feynman-diagram. For the calculation of 0(a 2 ) 
single bremsstrahlung corrections in this work, however, all physically outgoing 
particles will also be denoted by outgoing four-momenta, such that the scalar 
product of incoming with outgoing momenta will always be positive. Since 
massless spinors are used in this thesis, all terms homogeneous in m e will be 
dropped out of the expressions derived below. The masses, however, will be kept 
as regulators of "unphysical" collinear divergences inside the scalar integrals. 



2.2 Complex Loop Integrals and Bilinear Co- 
variants 

When using algebraic manipulation languages to calculate higher order Fcynman 
diagrams, it is important to reduce the sometimes very involved matrix expres- 
sions in numerators of the resulting integrals. Since these terms will eventually 
be evaluated between spinors, it is desirable to reduce strings of Dirac gamma 
matrices — up to seven in integrals that will be treated in chapter [| — to terms 
multiplying bilinear covariants . The transformation behaviour of the various 
basis terms is that of a scalar (1), a pseudoscalar (75), a vector (7^), an axial- 
vector (757^) and that of a second rank tensor {0^ = |[7^,7iy])- Expanding 
an arbitrary string of slashed four-vectors, G =^\... $ n , yields: 
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Figure 2.1: The t-channel tree level Feynman diagrams contributing to the 
0(a) single bremsstrahlung cross section for initial and final state electron line 
emission. 



G = \Tr{G)l + \Tr( l5 G) l5 + ±Tr{^Gyf 
-1^(757^)757^ + \Tr(a^G)a^ 



(2.36) 



It can be easily checked that if G is one of the stated basis vectors, Eq. [j.36j 
holds in each case. The use of 75 can be avoided in n-dimensions, but the alge- 
bra in the end always assumes that the spinor space remains four-dimensional. 
The reduction 2.36 also allows the application of the entire spinology that was 
described in the beginning of this chapter, since the terms in 2.36 are easily 
evaluated between spinors. Ref. |24[ contains a program written in FORM that 
produces the presented reduction in n-dimensions and is employed in most of 
the other routines used in this thesis. 



2.3 The Tree Level Single Bremsstrahlung Am- 
plitudes 



The 0(a) t-channel single bremsstrahlung electron- line emission amplitudes are 
shown in Fig. |2.l[ In the high-energy limit <C 1, one can legitimately make 
use of the massless spinor formulation described in chapter || and of the helicity 
conserving "magic" polarization- vector in Eq. 2.12] . The tree-level amplitude, 



from now on denoted by A T e - , can, for initial state electron- line emission, Fig. 



2.1, be expressed as 
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Al L 



i e 



:<Q',A'| 7/ ,|P',A'>Ga,*«) x 



-2PK 

(Q,A|y(p-jf)^|P,A) 



A,-P 



(K,P) X (K,P)_ X 
+ 2(Q',A / |^|P / ,A')(Q,A|^|P ! A) 

ie 3 y^pG A ,y(t / )^,- P , 
(K,P) x (K,P)_ x (h,K) p X 

{Q,P). X (P,K) X 



2Ps{Q',X'\^\P',X'){Q,X\^\P,X) 



(2.37) 



(2.38) 



(Q',Q)^ X ,{P,P') X , ,X 
{Q',P}.y{Q,P') x , ,X 



X' 



-A' 



+ (Q,^}_ A (X,P) A 

ie 3 y / 8flG A , A ft / )(5 A ,_ p 
(P,K)- X (h,K) p 

(Q,P)-x 



(Q',Q)- X ,(K,P>) X , ,A 
(Q',K)- X ,(Q,P') X , ,X 



X' 



-A' 



(Q',Q)- X ,(P,P') X , , A = A' 
(Q',P)- X >(Q,P') X , ,X = -X' 



(Q',Q)- X ,(K,pi) x , , A = A' 
(Q',K)_ X ,{Q,P') X , ,A = -A' 



- (Q,i^)-A 

»e 3 V8 P G A:A (^)(5 A ,_ p 
(P,K) p (h,K) p 

(Q>, Q) P ((Q, P)-x{P, P')x - (Q, K)-x{K, P') x ) , A = A' 
(Q, P%«Q, 0-a(Q', P)a - (0, #>-a(Q', *Qa) , A = -A' 

ieVSpG^^O^p UQ',Q) P {Q,Q')- X {Q',P')x ,A = A' 



(P,K) p (h,K) p \(Q i n p (Q 1 P')- X (P',Q , ) x ,A = -A' 



{h,K) p {h,K) p 



(2.39) 



where the propert ies of of Eq. [2.12 , four-momentum conservation and the 
Fierz-identity Eq. 2.9 were used. Eq. 2.39 , derived for a specific electron-photon 
helicity correlation, holds in fact for both hclicity cases. This was achi eved 
by using the helicity dependent /i-spinors of chapter ||. Although Eq. 2.39| 
represents the numerically most stable form for A T e - , there is another way of 
writing 2.37 by means of the identity 2.23] which will prove to be very useful in 
the calculation at a later point and is indeed the main application of 2.23: 
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Al L - = 



le 

hK 
2 



pXt - 2hK 



(Q , ,A , |7 M |P , ,A')<g,A|y|P ) A) 



(Q',\'W\P',\')(Q,\\fi\P,\)}he 



(2.40) 



Eq. 2.40 will play an important role in separating out the "true" tree-level part 
of the complete 0(a 2 ) internal emission calculation in chapter see also Ref. 



Chapter 3 



The 0(a 2 ) Single 
Bremsstrahlung Amplitudes 



This chapter is devoted to describing in detail the steps of the second order cal- 
culation. While a lot of analytical work is presented, all occuring loop integrals 
were evaluated in the appendices of Ref. Q with the algebraic manipulation 
language FORM. It is assumed from this point on that all UV-divergences will 
be properly taken care of by the appropriate counter terms in the Lagrangian 
such that in the on-shell renormalization scheme — e and m e correspond to the 
physically observable charge and mass of the electron respectively. 
Employing standard Feynman techniques |2l], [l2| gives 18 electron line emission 
graphs for the desired level of precision in the previously described context 
(compare with chapter |l|) . The single bremsstrahlung amplitudes contributing 
to 0(a 2 ) radiative corrections to the low angle Bhabha cross section are listed 



in Fig. 3.1 for t-channel electron-line emission. Positron-line emission can be 



obtained from these amplitudes through crossing relations: 

P <— > -Q' (3.1) 
P' < — > -Q (3.2) 
A — ► - A (3.3) 

A crossing rule can als o be used to calculate the initial state s-channel contri- 
bution by applying |3.2| to the electron-line t-channel calculation. The s-channel 



final state contribution is attained by employing 3.1 and 3.3 to the expres 



sions corresponding to the diagrams shown in Fig. 3.1. All vacuum polarization 
graphs that would in principle also belong to this class of diagrams are explicitly 
omitted here, since they are already included in BHLUMI4.01 |35||. 
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18CHAPTER 3. THE 0(a 2 ) SINGLE BREMSSTRAHLUNG AMPLITUDES 




Figure 3.1: The one particle t-channel exchange Feynman diagrams, modulo 
vacuum polarization graphs, for initial and final state electron line emission that 
give a nonzero contribution to the 0(a 2 ) single bremsstrahlung cross section 
after rcnormalization is carried out in the chosen on-shell scheme. Also Z- 
exchange can be included. 
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Following the notation of Fig. 3.1 the various graphs translate into the ampli- 



tudes given by the expressions listed below. It should be noted that the Feyn- 



man gauge [|23|, 42 was picked to arrive at the following representation. The 
divergent integrals are regularized by means of the n-dimensional regularization 
technique. 

^d n i (q, a] jjj t \p, a) 

-2PK (I 2 - m n 2 )((Z + P-K) 2 - m e 2 )((l + Q) 2 - m 2 ) + ie 

(3.4) 

M = -^(Q',\'\^\P',\')G x , x ,(t')x 

^d n i (q, Aj j @ +ph„u +0 + w(/ +fb ,y \P, A) 

2QK (I 2 - m 2 )((l + P) 2 - m 2 ){{l + Q + K) 2 - m 2 ) + ie 

(3.5) 

A 3 = ^(Q',A'| 7 jP',A')G A ,*(0x 

d^(Q ; A|7,(/+^)7 M (/+f-^)^/+f)7' y |^A) 
(I 2 - m 2 ){(l + P) 2 - m 2 ){{l + P-K) 2 - m e 2 )((l + Q) 2 - m 2 ) + ie 

(3.6) 

A, = -^{Q 1 ,\'\ ltl \P' ,\')G x , x {t') x 
( 27T ) 

d 4 i(Q, aj yjj+m iu+g j^ru+ryf \p, a) 

(l 2 - m a 2 ){{l + P) 2 - m e 2 ){{l + Q + K) 2 - m e 2 )((l + Q) 2 - m 2 ) + ie 

(3.7) 

(2tt) 4 



-2PK J ^ (I 2 - m 2 )((l + P - K) 2 - m e 2 ) + ie -2PK fl ' ' 

(3.8) 



A ^ = 7^(Q'^'h„\P',X')G x , x >(t')x 



(27T) 

2QK J 1 (I 2 - m 2 )((l + Q + K) 2 - m 2 ) + ie 2QK 1 ' ' 

(3.9) 



e 5 



A ? - -^{Q 1 ,X'\^\P' ,X')G X , X (t')x 



(2tt)' 

f-ff f fJ^ljuU+l -P)tU+f)l v |p \\ 

-2PK J {I 2 - m 2 )((l + P) 2 - m e 2 )((l + P - K) 2 - vn 2 ) + ie ' ; 

(3.10) 
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e 5 

(27T) 



^8 = 7^4<Q'> A 'l7^',A')GA,A'(t')x 



W ' 7 (Z2_ mo 2)(( / + Q)2_ me 2)((z + Q + x) 2_ me 2) +ie 2Qif 7 1 ' ; 

(3.11) 

^9 = ^!(Q,A| 7Al t|^P,A)G A ^(0x 

(z 2 - m 2 )((/ + ^') 2 - m e 2 )((; + go 2 - m e 2 ) + it 

(3.12) 

^d n i(Q\ \'\ lv (/+$')r(/+r)r\p', x) 

(I 2 - m 2 )((/ + P') 2 " m e 2 )({l + Q') 2 - m 2 ) + le 

(3.13) 

The following sections of this chapter will now go about solving the integrals 



appearing in Eqs. — 3.13 and give concise expressions for the corresponding 
amplitudes. 



3.1 Vertex Corrections 

The name "vertex corrections" to the 0(a) single bremsstrahlung cross section 
implies that some kind of a factorization takes place, that separates the lower 
order amplitudes from the additional virtual corrections. This, however, is not 
true in general and only holds in this calculation for A±, A2, Ag and A\q because 
of the application o f the high energy approximation and the use of th e "m agic" 
polarization vector [2.15 and the "magic" choices given in Eqs. |2.13| and Ell 



The graphs calculated in this section include amplitudes 3.4, |3.5|, 3.10, 3.11 



3.12 and 3.13. For the amplitudes A-j and A$,, the described factorization does 
not take place despite the tools mentioned above. It still will be useful in this 
case to use the broad headline "vertex correction" because of two reasons: 
The result for A7 and As is completely determined by form factors derived in 
this section and secondly, the terms not proportional to the lower order, tree 
level, amplitude turn out to be canceled exactly by a contribution contained in 



the rather complicated expression for A3 and A4, to be discussed in section 3.3 
This circumstance will insure that all leading and subleading contributions will 
be concentrated in the factorized part of the complete answer. 
All results employ the on-shell renormalization scheme, i.e. they are renormal- 
ized such that the vertex correction contribution vanishes if both fermion legs 
are on the mass-shell and for zero four-momentum transfer. 



3.1. VERTEX CORRECTIONS 



21 



3.1.1 Scalar Integrals 

The Scalar integrals, in terms of which the following results of the various am- 
plitudes are given, are defined in this section as follows: 



£>12 
B23 
£>23 
123 
Cl230 



B (m e 2 - a;m ,m e ) 

B (2m^+t'-m e ,m e ) 

B (m 2 ;m e ,m e ) 

Co(m 2 - a, 2m 2 +t', m 2 ; m ,m e ,m e ) 
C (trie , m Q 2 , m e 2 , m , m e , m e ) 



(3.14) 
(3.15) 
(3.16) 
(3.17) 
(3.18) 



It must be emphasized that the notation for scalar integrals in 3.14 — 3. IS only 
refers to loop integrals with no more than three internal propagators. This 
was done in order to have agreement with the notation implicitly employed by 
FORM in the appendices of Ref. @. 

Eq. [3.18j is somewhat special among the integrals denoted above in that it can 
be expressed analytically and without renormalization and scale dependence: 



C 



1230 



i I <I'I 



1 



(I 2 ~m 2 )(l 2 + 21P) 2 + ie 



dx- 



x 



it:- .1 J (I 2 + 2lPx - m 2 (l - x) + ie) 3 
• 1 

dx x 

11 m 2 x 2 — mtfx + Too 2 + ie 



2m,_ 
vnr 2 In 



1 1 

7 / 22 2 1 2 1 • \ 
ln(m e x — toq x + mo + ie) 



toq 



(3.19) 



where standard Feynman parameters |42j, |2l| , integral tables J|14| and formulas 
for integrals over four- dimen sional Minkowski space |2l| , [27| were used. In 
addition it follows from [3.18 that P 2 = m 2 : so that only masses are parameters 
of this integral. The significance of [3.19| here is merely a notational one. It 
means that m e 2 Ci23o-terms are actually of order O(m e ) and need to be taken 
into account; it also is IR-divergent, which implies the possibility of a strong 
analytical check of the expected limit for mo — > 0; see chapter]^ for an analysis of 
this problem. The same behaviour with respect to a possible m e ~ 2 dependence 
is not given for the other integrals listed. Dimensional reasons exclude it for 
Aq and Bo functions and numerical checks for the other integrals in the regime 
where products of four- momenta are much larger than m 2 . 



3.1.2 The Vertex Function 

The vertex function Tg((P — Q) 2 , P 2 , Q 2 ) is given by the integral: 
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H e d n l 



(I 2 - m 2 )((l + P) 2 - m 2 )((l + Q) 2 - m 2 ) + le 



(3.20) 



The numerator of the integral 3.20 can be simplified by means of standard 
Dirac-algebra manipulations as follows: 



= 7,(/ W + m e 7„7"(/+ W + m e7 „(/+$)yV 

= -2(1 +f ) 7 "(/ +<$) + 4m e (l + PY + 4m e (i + Q)" - 2m* -y" 

= -2((- 7 "(/+P) + 2(1 + Pf)(l+Q) + (m e - terms) 

= -4(/+#)0 + P)" + 2 7 "(-(/ + $)(/+P ) + 2(1 + P)(l + Q)) 

+ (m e — terms) 
= -4(/+0)(/ + P)» + 4rf(l + P)(l + Q) + 2((/+^) 7 " 

-2(/ + Qy )(/+]/>) + (m e - terms) 

= 2(/ (; +?) + 4 7 " (j + p)(i + q) - 4(i + pyu+Q) 

-4(1 + QY l (l+f ) + (m e - terms) 
= 2{f+Q - m e )f(/ +P - m e ) + 2m e7 ^(/+f ) + 2m e (/+#) 7 /i 

+4rf{l +P)(l + Q)- 4(1 + Py(f+Q) - 4(1 + Qy(/+f) 

+4m e (l + Py + 4m e (l + Qy - 4m e V 
= 2(/+0 - m e )^(l+P - m e ) + 2m e7 ^(/+f - m e ) 

+2m e (/ + $ - m e )y + 4^(1 + P)(l + Q) - 4(1 + Py(/+@- m e ) 

-4(l + QY l (/+P-m e ) 
= 4^(1 + P)(l + Q) + 4m„r - 4 /(2i + P + QY - 2^l 2 

+4l"(/+@ - m e ) - 2(0 - m e ) / 7 " + 2 / 7 "(P - m e ) 

+2(0 - m e )y(P - m e ) + 2m e7 /i (P - m e ) + 2m e ($ - m e ) 7 " 

-4(1 + PY(Q - m e ) - 4(1 + QY(P - m e ) 
I 2 

= 4^(PQ + l(P + Q) + —)-4/(P + Q + iy + 4m e P 

+ (<# - m e )(2m e Y -IJrf- 4(1 + P)^ + 41") 

+ (2m e7 " + 2 / 7 M - 4(1 + QY)(f - m e ) 
+2(Q - m e ) 7 "(f - m e ) 

- 4 7 ^(PQ + l(P + Q) + l -) - 4 /(P + Q + iy + 4m e i" 
+ (<# - m e )(2m e7 ^ + 2 7 " / - 4(/ + P)^ 1 ) 
+ (2m e7 M + 2 /7" - 4(1 + QY)(P - m e ) 

+2($-m e ) 7 M (P -me) (3.21) 
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The first line of Eq. 3.21 contributes regardless of whether the fermion lines 
are on- or off-shell and can be found in standard textbooks, for example in Rcf. 
p9[. In case one fermion line is on- the other one off-shell, either the second 



or third line, respectively, contributes in addition. The last term in Eq. 3.21 



is only needed for the case where both P and Q are off-shell and thus will not 
enter into the present 0{a 2 ) level calculation. The mass-terms are retained only 
until the last stage of the FORM-reduction, after which only terms of O(m e ) 
are kept; see the //ac-routines in Ref. [E4J for details. 



For initial state electron-line radiation, Eqs. 3.4, 3.12, the renormalized vertex 
function V 1 can be written as 



r n = FFY L + FF a {f -fi)(P - KY + FF b (tf -ft)Q^ 



(3.22) 



With the previously defined notation for scalar integrals in section 3.1.1 for the 
external emission graphs A\ and Ag, the decomposition 3.22 is done by the 
FORM-routine ffac, Ref. p4| , and yields the following results for the various 
form- factor terms up to O(m e ): 



FF — —2 — 4m e 2 Ci23o — 2i'Ci23 + B 23 + 2Bf 2 — 3i?23 
2 

FF n = (B-\o — 

FF b = -4C? 23 +-^(B^ + l + aCf 23 ) 
t + a I 

+ , 2 (t'B?2- \aBf 2 +\aB 23 - 2t'B 23 ) , 



(3.23) 
(3.24) 

(3.25) 



for the off-shell graph A\ and 



FF° 

FFa 



FF?, 



2 - 4m e 2 Ci23o 


0. 



2i'Ci23 + 3i?°3 — &B 



>2.\ 



(3.26) 
(3.27) 
(3.28) 



for the on-shell graph A 9 . 

The identity (f -$T) = + 2{f -fi){P-Ky was used in the above 
decomposition as well as the renormalization condition that vanish on-shell 
and for zero transfer, t' = 0. Due to the choice of the polarization-vector e M 
the expression (Q,A\$\P, A) = 0, which means that the form-factors FF a and 
FF b don't contribute to the aspired level of accuracy for Ai{A 2 ) and Ag^io). 
They will, however, prove to be very useful for te sting the internal consistency 
of the calc ulation, as will be discussed in chapter Qf2 and do contribute to the 
amplitudes 3.1C| , |3.11 . While it is not directly apparent, it can easily be checked 
that all renormalized form factors are indeed finite and scale independent. Ref. 
\QM also demonstrates that the renormalization condition is met correctly. 



2ACHAPTER 3. THE 0{a 2 ) SINGLE BREMSSTRAHLUNG AMPLITUDES 



Inserting Eq. 3.20 into the expression for A\, Eq. 3.4, and using the decompo- 
sition Eq. 3.22 as well as the "magic" properties of e M yields: 



A l = 



(Q',X'\^\P',X')G x , x (t') x 



16tt 2 

(Q,A|r"((P' -Q') 2 ,(P-K) 2 ,m e 2 ) 



-2PK 



\P,X) 



ie 



16tt 2 

{(Q,Air 



(Q',\'\^\P' : \')G x , x ,(t') x 

JP - Tt 

P, A) FF 



-2PK 

(Q, A|(f -$f^\P, A) (FF a (P - KY + FF b Q») } 



ie 

16n- 



{Q',X'\^\P f ) X')G x ,>i(l/) x 



(QA\rz^$\PA)FF 

A T e - FF{{P' — Q') 2 , (P — K) 2 ) 6 p> -\ 



16tH 



For Ai one gets analogously: 

„2 



A, = 



16tt 



2 A T e hFF{{P' -Q') 2 ,(Q + K) 2 )S, 



(3.29) 



(3.30) 



Since the tree-level amplitude A T e - already contains the helicity dependence 
implicitly it is convenient to write: 



with 



A x + A 2 = 



7 



16tt 2 



A T ^FF(t',^, 



m e 2 + 2QK , p = A 



2PK 



-X 



(3.31) 



(3.32) 



For the on-shell form-factor amplitudes Ag and Aio the derivation is even more 
straightforward and gives: 



A 9 + A w = 



167T' 



■A T e ±FF°{t') 



(3.33) 



As can be seen from the definition of the scalar integrals for the various form- 
factors the relation between FF and FF° is given by: 



FF°(t') = FF(t'^ = m e 2 ) 
This correspondence is shown to hold numerically in Fig. |3.2|. 



(3.34) 
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Form Factor Comparison 




Figure 3.2: The real parts of the on-shell (dashed lines) and off-shell form factors 
for various values of -^-r — 1 at t' = —1. The off shell form factor is seen to 

me 

approach FF° for 7 — m, c ? 



0, Eq. 3.32, and to have no log— divergence for 



7 > m e 



3.1.3 The Self-Energy Emission Amplitudes 

The contribution from the off-shell self-energy emission amplitudes 



can be described by the expressions for the off-shell form factors 3.23| , |3.24| and 
3.25 in the limit of zero transfer (t 1 = 0) as follows: 



3TTo| pd 



e^(K 2 ,P 2 ,{P-K) 2 ) = 



FFi + {FF a 
(-2 - 4m e 2 C 1230 + 



23 ~ -B12 ) i + 



•4B§ 3 -4Bf 2 )(p-^) (3.35) 



Using the above limit and decomposition gives for amplitude 3.10 



A- 



16tt 2 

5 



-2 - 4m e 2 C 1230 + 5 2 ° 3 - Bh)A T e h 



IP" F PL 1 

-^(6 + 4B 2 °3 - 4Bf 2 )^- G x , M (t')(Q', A' I 7p \P', A') (Q, A| 7^ \P, A) 
1D7T a 
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16tt 2 

=5 



(-2 - 4m e 2 C 1230 + B 2 °3 - B?2 ) A 



TL 



ie 



(6 + 4B 2 ° 3 -4B? 2 ) 



t + 



(Q',\'\Q\P',\')(Q,\\fi\P 7 \) - 
t-0 



Gx,x(t') x 



atA 



TL 



16tt 2 



ie 



t + 



2e^ S 



[Bh - S 2 ° 3 ] - 4m e 2 Ci230 



2ie 3 (t + /3) 

5f + 2/3 



X-p 



16tt 2 (t + /3) a 
(Q',\'\Q\P',\') (Q,\\£\P,\) 



(6 + iB$ 3 -iB? 2 )G x ,x(t') x 



.4 



TL 



(3.36) 



The important identity 2.40 was used in deriving the result in 3.36, in which 
all terms proportional to A F e - were identified. As will be discussed in more 
detail in chapter 7.3, terms multiplying the 0(a) single bremsstrahlung tree 
level amplitude contain the biggest corrections to the 0{a 2 ) cross section. The 



result for A s can be obtained from Eq. 3.36 through the crossing rule 2.22 



A 7 (P^ -Q,X 



-A) 



(3.37) 



3.2 The Electron Self-Energy Function 

|, are determined by 



The contribution of the amplitudes given by Eqs 
the solution of the electron self-energy function 



£ o (P,TO ,TO e ) = [ n e d n l- 

Z7T 2 J ( 



(3.38) 



{I 2 - m 2 ){{P + I) 2 - m e 2 ) + ie 
K(P 2 ,m ,m e )m e + £ b (P 2 , m , m e )f> (3.39) 



The proper treatment of the mass terms in the expression 3.3S for the self- 
energy contribution is essential here, since there are two independent renor- 
malization constants in £ (P, mo, m e ). In keeping with the hitherto adopted 
on-shell renormalization scheme, these two constants are fixed by imposing the 
renormalization conditions 



£(P, m ,m e ) 
<9X(P, TO ,m e ) 



dP u 



= 



(3.40) 
(3.41) 



The renormalization condition |3.41 is the reason why the terms proportional 



to m K matter here: §§-L_ contains terms proportional to — . After prop- 



3.2. THE ELECTRON SELF-ENERGY FUNCTION 



27 



erly taking these into account, the FORM-routine selfen, 
following result: 



Ref. 



gives the 



£(P, to , m e ) = [4 (B? 2 - P? 2 ) - 4TO e 2 C 1230 - 2] to 
4TO e 2 C 1230 - (P? 2 - B° 12 ) M + ^ 
» [4TO e 2 C 1230 - (Bf 2 - S? 2 ) + 3] f 



TO e 



(3.42) 



The mass-terms can be dropped in the final step in Eq. |3.42| since it represents 
the renormalized result. Doing the same at an earlier stage, however, would 
give the wrong S f, (P 2 , too, m e )-function. Eq. 3.42 is also clearly finite and scale 
independent. The "Co" function enters the self-energy result, a term with only 
two propagators, through the derivative in 3.41 and is given explicitly in Eq. 
3.19. Inserting this result into Eq. 3.8 yields: 



A 5 = 



16tt 2 



(Q',A'| 7ai |P',A')Ga,a'(0 x 



(Q, A| l^^Jr; S 6 ((P - if) 2 , too, to c )(P - Jf) I— % } |P, A) 



-2PK 



-2PK 



S 6 ((P - K f, TOO, TO e ) (Q', A'| 7m |P', A')Ga,*(0 x 



16?r 2 

(Q,A| 7 ^|^|P,A) 

„2 

E b ((P-if) 2 ,TOo,TO e )A 



i'L 



167T* 



(3.43) 



With the notation of Eq. 3.32 one finally gets: 



A 5 +Ae = 



167T 



S fc (7,TO ,TO e )A 



Ti 



(3.44) 



arc 1 



It is interesting to observe, that the IR-divergences contained in Eq. [3.44 
exactly canceled by the In^ 1 terms present in the expression for the amplitudes 
A7 and Ag in Eqs. 3.36 and 3.37 . This set of graphs is also separately gauge 
invariant. 

Standard Feynman techniques give 18 one photon exchange diagrams for the 
process considered here (modulo vacuum polarization) as was mentioned in sec- 
tion ||. From Eqs. 3.4C , 3.41 it can readily be seen that all graphs with self- 
energy loops next to external lines don't contribute to the 0(a 2 ) calculation in 
the chosen on-shell renormalization scheme, since 



S(P, TOo, TO e ) ^=m e <9E(P, TOo,TO e ) /=m, 
P-TO e ' Hp ' 



(3.45) 



Thus 
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A n = , ... , A 1S = (3.46) 

This completes the O(m e ) contributions of the 0(a 2 ) single bremsstrahhmg 
amplitudes with real hard photon emission off external legs. The next section 
will now deal with the much more complicated situation given for internal pho- 



ton emission in amplitudes Eqs. 3.6, 3.7 



3.3 The Internal Emission Amplitudes 
3.3.1 Scalar Integrals 

The FORM notation of scalar integrals used in the decomposition of the tensor 
integrals 3.6, 3.7 differs from the notation used in previous sections for the 
obvious reason that there are now four propagators inside the loop. This leads 
to the following identifications: 



Bl2 


= B a (m e 2 ;m ,m e ) 


(3.47) 


r>a 
£>13 


= B (m e 2 - a;m ,m e ) 


(3.48) 


#24 


= B Q (2m e 2 +t;m e ,m e ) 


(3.49) 


B 3 4 


= B a (2m e 2 +t / ;m e ,m e ) 


(3.50) 




= Co(m e 2 , m 2 , m e 2 - a; mo, m e , m e ) 


(3.51) 


C124 


= C (m e 2 ,2m e 2 +t,m e 2 ;m ,m e ,m e ) 


(3.52) 


W34 


= C (m e 2 - a,2m e 2 +t',m e 2 ;m ,m e ,m e ) 


(3.53) 


C234 


= Co(mo 2 , 2m e 2 + t' , 2m e 2 + t; m e , m e , m e ) 


(3.54) 


-Dl234 


= D (m 2 ,mo ,2m 2 + t',m 2 ,m 2 — a,2m 2 + t;m 07 m e ,m e ,m e ) 


,(3.55) 



where the basic B , Cq and D functions are defined in section 2.1 as usual. 
3.3.2 The A 3 + A 4 Calculation 

The nu mera tor of Eq. |3.6| was put into the FORM-f ile a3 a4 in which the outgo- 
ing term 3.7 was added through the crossing relation 2.22 . By means of algebraic 
manipulations, the initial output of roughly 20000 (!) terms could be reduced 
to 90, which involved many rathe r com plicated transformations (see appendices 
of Ref. (2^]). Again, the identity 2.23 was crucial to achieve this step down to 
so relatively few terms as well as separating out the part proportional to A T e - . 
The result given by a3a4 reads 



A 3 + At 



le 
16ir 



where 



; G x ,x>(t') (F 1 + F^i + F 2 1 2 ) , 



i = 2^,0-, 



(3.56) 



(3.57) 
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X, = W^^/f'^ (3.58) 
(h,K) p (P,Q)- P 

1 2 = 2V - 2X %J^^llf^l (3.59) 

(h,K) p (P,Q)- P 

The function To is proportional to the electron-line tree level amplitude: 

A T e - — ie 3 G\ t y(t') To (3.60) 

The form factors are 

F (p = A) = -2(«7i24 - t'Cf^) - (fl^ - B 34 )3/?(t - a)" 1 

{tC 124 - aCf 23 - (t + /?)<7?34 + (a - /?)C 234 + crt£>? 234 } 
x tp-^a- /3)(t- a)' 1 
+ {tC 12i + pC-£ - (t - a)Cf^ + (a - /?)C 234 - W2L I 
+ (B? 3 - B 34 ) a(t - a)- 1 {1 - 3t(t + p)- 1 } 

+ (B 24 - B 3A ) 2ta{a - [3)-\t - a)' 1 + (B 12 - B~f)2t(t - a)' 1 

(3.61) 

F 1 (p = X) = 2t(a-f3)- 1 

+ {tC 124 - aC? 23 - (t + p)C? 34 + (a - /3)C 234 + crt£>« 34 } 
x{tt'f3- 2 (t - a)- 1 ^ -(3) + \8p.x] 

- {tC 124 + 0Ci& - (t - a)Ci& + (a - /?)C 234 - /ftD^ } 

x{ia _1 /3^ P ,-i} 
+ (B~ 3 -B 12 ) 2t{t-a)- 1 
+ (Bf 3 - B M ) tt'(t - a)- 1 {2/T 1 - 3(t + !3)- 1 } 
+ 2(B 2A - B 34 ) tt'(a - ^{(a - /3)- 1 - tp-\t - a)' 1 } 

(3.62) 

F 2 (p = A) = -2t(a - (3)- 1 + t(t + (3)- 1 

- {tC 124 - aC? 23 - (t + /3)C? 34 + (a - p)C 23i + atD<? 23i } 

x{tt'l3- 2 + \ap- 1 5 p . l } 

+ {tC 12i + pCii - (t - a)Ci& +(a- /3)C 234 - ptD^ 3i } 
x{^p.-i} 

+ (S 34 - B? 3 ) tt'(t + P)- 1 {2/T 1 + (t + P)- 1 } 
+ 2{B 2A - B 34 ) tt'(a - (3)- 1 {/T 1 -{a- PY 1 } 

(3.63) 

The opposite helicity cases may be obtained from the above results using the 
substitutions (for i = 0, 1, 2): 
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Fi(p = -A, a, (3) =F l {p = A, -(3, -a) (3.64) 

While Eq. |3.56| still does not look very transparent, several observations may 
help to clarify the overall picture. Since the amplitudes |3.6| and |3.7| are not UV- 



divergent one should evidently expect this feature for the result 3.56, where all 
Bq terms possess such an infinite contribution (see section 2.1). A simple check 
in the FORM routine a3a4, which replaced each Bij by Bij +div, proved that all 
divergences cancel properly. This procedure also insures the scale independence 
of Eq. 3.56. Numerical analysis revealed that in the case of amplitude A3, for 
instance, only C124 and Z)"234have IR-divergent terms, regulated by the virtual 
photon mass mo. By closer examining the four denominators in 3.6, it becomes 
clear that , except for K° = 0, only propagators 1, 2 & 4 taken together and with 
no loop momentum left in the numerator can lead to a logarithmic divergence 
in mo, confirming the numerical findings. From Ref. |3l| and for r -^- -C 1 one 
has: 



C124 
and for mo — > 0: 



-lr 2 



-t 



+ In — 



wo 



7T~ 
"if 



(3.65) 



234 



at 



-In- 



-In 



mo 



(3.66) 



Eqs. 3.65j and 3.6(: corroborate that the only effective IR-divergent term in 



3.5E is — 2fCi24 iTpj-Zo) which coincides with the prediction of the IR-theory 



expounded on in chapter [|. It should also be noted that the form factors F\ 
and i*2 are normalized in a slightly different way than the output printed in 
Ref. p4| . For comparison, the form factors here should be multiplied by ^ for 

p = A. Pi contains a contribution proportional to (B 13 ^ — B1 2) w hich will turn 



out to cancel the corresponding contribution from amplitude 3.36| and will thus 
make F\ and F2 almost symmetric in form. This feature will help concentrate 
the numerically relevant leading and subleading logarithmic terms in the part 
proportional to the tree level amplitude, as the results of this chapter will now 
be summarized in chapter W. 



Chapter 4 



The Exact Differential 
0(a 2 ) Result 



In order to have a complete differential result it is necessary to express all 
amplitu des w ith the same notation of scalar int egral s . Using th e defi nitio ns of 
section 3.3.1 and the expressions given in Eqs. 3.31 , |3.33 , 3.36| , 3.37| and 3.44 
one gets the following for p = — A: 



Ai + A 2 + A 5 + A 6 + A 7 + A$ + A g + A 



10 



16tH 



P , R 4t + 6/3 j-y a 2i 3a 

-t>12— " r -Di 



t + P 



13 t' + a 



t' + a 



(4.1) 



and thus for p = A: 



„2 



An + A* 



16tt 2 
B34 



A 7 + A s + A g 
At - 6a 



B 



6*' ~ 3/3 o _ 2 ^ 



/3 



16tt 



t — a 

-0 



t — a 



D 



A w 



13 



t'-P 



zr u 134 o 134 



2B 



13 



2B, 



I 



(4.2) 



The exact C(a 2 ) single bremsstrahlung result for electron -line emission and 
p = A is of co urse given by adding all amplitudes 3.4 + ... + 3.13 and with Eqs. 
3.56 and \L2 reads for p = A: 
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A P - = 



16n' 



;G X ,\>{t') (Wo + Till + Till) 



(4.3) 



where the notation of Eqs. 3.57, 3.58 and 3.59| has been used. The form factors 



arc 1 



T (p = A) = -8 + a(t - a)- 1 - 8m e 2 C? 23 - 2(tC 124 + i'C? 34 ) 

- {tC 124 - aC? 23 - (t + [3)C? 34 + (a - f3)C 234 + atDf 234 } 

x tfi- 1 ^- P)(t- a)' 1 
+ {tC l24 + f3Cii - (t - a)Ci£ +(a~ {3)C 234 ~ ptD^ } 
+ 6B 12 + (B° 3 - B 34 ) a(t - a)- 1 {1 - M(t + /3)- 1 } 

- 6B 34 + (B 24 -B 34 ) {2ta{a- f3)- 1 {t-a)- 1 } 



Ti(p = A) = 2t(a - - t(t - a)' 1 

+ {tC 124 - aC? 23 - (t + f3)C? 34 +(a- f3)C 234 + atD« 234 } 

x{tt'f3- 2 (t - a)-\t -/3) + i<5 p ,i} 
- {tC 124 + pC<i£ - (t - a)Ci£ + (a - fi)C 234 - PtDiL } 

xjia- 1 /?^-!} 
+ (Bf 3 - B 34 ) tt'(t - a)- 1 {2/3- 1 - 3(t + /3)- 1 } 
+ 2{B 24 - B 34 ) tt'(a - ^{{a - /3)" 1 - tpr\t - a)' 1 } 



MP = A ) = -2*(a - PT 1 + t(t + (3)- 1 

{tC 124 - aC? 23 - (t + /3)C? 34 + (a - (3)C 234 + atD? 234 } 
x{tt'f3- 2 + \a[3' 1 5 pA } 

+ {tCi 24 + PC^ 3 - (t - a)C u 1 +(a- f3)C 234 - ptD^ 34 } 
x{^p,-i} 

+ (B 34 - B« 3 ) tt'(t + /3)- 1 {2/T 1 + (t + /3)- 1 } 
+ 2{B 2A - B 34 ) tt'(a - P)' 1 {/T 1 - (a - PY 1 } 



(4.4) 



(4.5) 



(4.6) 



The opposite helicity cases may again be obtained from the above results using 
the substitutions (for z = 0,1,2): 



Ti(p = -A, a, (3) = Ti(p = A, -0, -a) 



(4.7) 



In Figs. 4J and 4.2 the differential result |13J is compared to a recent calculation 
by Fadin et al. [|43f| after the IR-divergent terms were properly removed (see 
chapter ||) . 
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virtual photon program 



July 24, 1995 



Run 01 



Beam energy — 45.6 GeV 
Kinematic mass corrections included 
Trigger thresh hold — .0xs B 
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Figure 4.1: A comparison of the exact differential expression |4.3| to the re- 
sult given by Fadin et al. p3[. The graph on top displays \M.°( a+a )| 2 — 
2a[Re(B{t))+Re(B(t , ))]\M°<^>\ 2 (see chapter |) , i.e. there is no IR-dependence 
(mo or E cu t) left over. M. denotes the invariant matrix elements for the respec- 
tive cases. The second figure shows the ratio of the two calculations explicitly 
over the range of available photon energies and the bottom graph gives the size of 
the 0(a 2 ) virtual corrections to the differential 0(a) single hard bremsstrahlung 
cross section. 
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virtual photon program 



July 24, 19 



Run 07 



Beam energy — 45.6 GeV 
Kinematic mass corrections included 
Trigger threshhold — .OxE B 
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Figure 4.2: This figure is analogous to Fig 4.1 but shows the same compar- 
ison for a different set of photon and positron angles. While runl has the 
bremsstrahlung photon close to the incoming electron direction, this graph has 
the emitted photon along the incoming positron line. Although the respective 
lines are dominating in each case, both figures are using the contribution from 
the complete two line t-channel expression. 
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The agreement between these two differential results is roughly within 5% over 
the range of available real photon energies. This should be compared to the 
overall size of the second order correction to the Born process, which is about 
3.5pm Jf5| . The soft limit, however, disagrees, which will become more im- 
portant for the integrated cross section since it is the regime with the highest 
emission probability [^|. As can also be seen in both figures, there is no virtual 
photon-mass dependence left. The one graph that is shown is actually three 
identical graphs for three different values of too (compare with figure legend) . 
The general treatment of the expected IR-behaviour is discussed in chapter || 
in the context of Yennie, Frautschi and Suura theory [^J. The expression is ex- 
plicitly scale independent, as is expected of every physically meaningful result. 
This holds not only for 4.2 but also for all equations that have contributed to 
the final answer. It is also worth mentioning that the by far largest part of the 
answer 4.3 is concentrated in only 18 terms, namely those contained in J-q. Nu- 



merical tests showed that on the scale of Figs. 4.1 and 4.2 no noticeable change 
could be observed. Since these represent only one dimensional cuts through a 
multidimensional phase space, this aspect needs to be revisited in the context 
of Monte Carlo (MC) phase space integration. 

Eq. [D] is the central result of the presented work; it describes the fully differ- 
ential 0(a 2 ) single bremsstrahlung corrections to low angle Bhabha scattering, 
which is essential for any MC-implementation pS4| . It is thus of primary im- 
portance for the theoretical luminosity determination below the 1pm precision 
level. In order to convey a more in depth understanding about the validity of 
the presented work, part II of this thesis will mostly be devoted to demonstrate 
that Eq. 4.3 satisfies all theoretical expectations and is in good agreement with 
leading log calculations. Also Monte Carlo results for the fully integrated cross 
section will be presented for a SICAL-type acceptance [Eq, km . 
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Chapter 5 

Yennie, Prautschi & Suura 
Theory 



Although quantum field theories with massless gauge bosons display both short- 
(UV) and long-distance (IR) divergences in certain types of graphs, it has long 
been known that the nature of the two infinities have very different origins 
p3| , |l5| . While the UV-divergences are a manifestation of our ignorance about 
small-scale physics and cannot be avoided in theories that assume point-like 
elementary particles, the genuine infinite wavelength divergences should never 
be "real", i.e. should cancel out of the final result for any cross section. The 
cancelation of the first order IR-divergences in QED was first proven by Bloch 
and Nordsieck || as early as 1937. It took, however, almost another quarter of 
a century until the IR-problem was understood comprehensively to all orders 
in the fundamental work by Yennie, Frautschi and Suura (YFS) Q. Since the 
YFS-theory gives a well defined IR-limit for the amplitude 3.56| , it is necessary 
at this point to give a brief review of the IR-theory. a in this chapter denotes 
the fine-structure constant. 



5.1 The YFS B & B Functions 

Consider the expansion of the full connected amplitude for s-channcl Bhabha 
scattering at y/s — Mz in terms of the number of virtual photon loops. One 
can then write Q 

oo 

M{P,P') = ^M„(P,P') , (5.1) 

n=0 

where M n (P, P') is the contribution of all n virtual 7-loop graphs to M(P, P'). 
The result obtained by Yennie, Frautschi and Suura is that 
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M n (P,P') = J2 m n-r(P,P' 



(5.2) 



where the rrij(P,P') do not have virtual IR-divergences and are of order a? 
relative to Mq — ttiq- The virtual infrared function B(s) is such that 



B(s) 



d 4 k 



(27r) 3 J k 2 — tuq 2 + it 



-2p; - k M 



^Pfl kp 



-2P'k -k 2 -it 2Pk - k 2 - ie 



ReB(s) = -1 

Z7T 



Hence 



s / , m e 1, s l\ „, m e „ 27T 
/n — ;? ( 2/n— + -in — ^ 2ln— + 1 



mo 2 m e : 



m 



(5.3) 



M(P,P') = e aB ^m„(P,P') 



(5.4) 



This is the famous exponentiation of virtual infrared divergences of the YFS 
program. To complete the review of the YFS-theory, consider next the differ- 
ential cross section for the process e + + e~ — > 71(7) + X , where 71(7) represents 
the emission of n real photons with four- momenta k n . For a given value 

of n, this differential cross section is 



da = 



VaReB 1 f TT dE x d 3 kj 



n 



E 

7=0 



771; 



.(«) 



(5-5) 

where m/™^ is a special case of rrij in [5.4| in which X involves n real photons. 
The second theorem of Yennie, Frautschi and Suura is that 



1=0 



S(ki)...S(k n )Po + J2 S(k 1 )...S(k i - 1 )S(k i+1 )...S(k n )0 1 (k i ) + ... 

1=1 

n 

+ y)S(ki)/3 n -l(ki, ... ,fcj_i,fej + i, — ,fen) +/3„(fci, ...,/c„) , (5.6) 



where /3j is IR-divergence free and is of order a J relative to (3q- The real infrared 
divergence function S(k) is given by 



S(k) 



a 
4^2 



P P' 
kP kP' 



(5.7) 



It follows that the cross section for the emission of an arbitrary number of real 
photons can be written as 



5.1. THE YFS B &B FUNCTIONS 
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do = e 2a ( ReB+B ^— / dye l ^- Ex)+D 



71=1 J j = l 3 



dE 



x 



(5. 



with D and B defined as follows: 



D 



and 



k°<E cut rf3 . 



(5.9) 



B(s,E c 



k"<E c 



1 

2^ 



d 3 k S(k) 



m - 



, 2 2a 



lu- 



ll) r 



m e 1 s , £■£?' 

2m h -m — =■ — in—7r 

m 2 to/ ESat 



-2ln - 

m E£ ut 3 



(5.10) 



From the explicit forms in Eqs. and |5.10| it can be verified that Re B + B 
is free of IR-divergences so that da in 5^ is indeed a physically meaningful 
quantity and exhibits the cancelation of infrared divergences to all orders in a. 
It should be mentioned that the infrared functions given above are exp ress ed in 
the s-channel and that the t-channel expressions are derived from |5.3| and 5.10| 
with the analytical continuation technique and read |3j| : 



ReB(t) =-i 

ZTT 
1 

B{t,E cut ) = — 



In — 7T 2ln h —In — - 

z » mo 2 m e 



III,: 



2ln— + 1 - — 
Too o 



(5.11) 



i — - I 2m 1 — In — r — In — _„ I 

m 



In — k 2ln 1 — In — ~ — In- , , 

m£ V mo 2 m/ E£ ut J 



„, m e , EpEn 
-2ln— + In J 4 

m E c z ut 



7T J 



(5.12) 



It is with the r eal Y FS-IR function 5.12] that the infrared limit of the internal 
emission result 3.56| was tested. Numerical results for two representative kine- 
matical cases are shown in Figs. 5.1. As is explained in the figure captions, the 
result for the amplitudes A\ — A\q is shown to have the right IR-terms predicted 
by the YFS program. 
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Figure 5.1: For the two separate runs, specified in Figs. 11, 12 as l(R-l) 



and 7(R7), the ratios of the squared invariant matrix elements resulting from 

the exact differential expression p~3| versus 2a(B(t, E cut ) + B(t', E cu t))\M 0< - a ^\ 2 

are plotted. The limit mo — > should be k — E ^ K independent and the 

ratios shown should approach —1 since the real YFS-functions are used. Both 

constraints are shown to hold for the result O. The finite gap to —1 seen in 

l 



both cases is of order 



log(m ) 



Chapter 6 

Gauge Variations 



Ever since t'Hooft was able to prove the renormalizability of spontaneously 
broken gauge field theories |l6[ , gauge invariance has become the dominant 
principle in modern particle physics. While it was long known that QED has 
gauge symmetry [|[ |3(J, the weak and also the strong (QCD) nuclear force 
are now formulated through gauge invariance principles p3| , |4^| . Many of the 
consequences of gauge symmetries such as the Ward-Takahashi identities can 
be used as cross checks of results obtained from related amplitudes. One such 
cross check in the present calculation is a result of the gauge invariance of a 
subset of amplitudes and will be expounded on in the following sections. 



6.1 Gauge Invariance 

All physically observable quantities in the Standard Model must naturally be 
gauge invariant, however, also unobservable parts of those quantities can already 
display this symmetry. It can be shown that the amplitudes 
{A 5 + A 6 + Aj + A s } and {A9 + A10} are all separately gauge invariant. This will 
now be shown explicitly for the first bracket of amplitudes. Gauge invariance 
demands that longitudinal photons decouple from the amplitudes ]3C|] , thus by 
replacing j. — >]jC one gets: 

A l r a = -^^A'| 7 jP',A')G A ^')x 
( 27T ) 

^d n l (Q, A| 7„(; -JTrf (g -£)jr |P, A) 

-2PK (I 2 - 77i 2 )((7 + P-K) 2 - m e 2 )({l + Q) 2 - m 2 ) + ie 

(6.1) 

^d n i (q, A|^r {Q +^yu/+Q +$)Y{l+f)i v \p, a) 

2QK (I 2 - m 2 ){{l + P) 2 - m 2 ){{l + Q + K) 2 - m 2 ) + ie 

(6.2) 
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A l r 9 = j^(Q',X'h,\P',X')G x , x (t')x 

dH(Q, a; jjj +oyni +i -$W (/ +1 )l v \P, A) 

(Z 2 - m 2 W + Pf - m e 2 )((l + P - Kf - m e 2 )((l + Q) 2 - m e 2 ) + ^e 

(6.3) 

^4°" s = j^y 4 (Q',X'\i,\P',X')G^(t')x 

dH(Q, Aj jjj+m u+q \p, x) 

m 2 )((l + P) 2 - m e 2 )((/ + g + X) 2 - m e 2 )((Z + Q) 2 - m e 2 ) + le 

(6.4) 



Using the identities 

u+r-pwu+p) = (i+p) 2 (/+p-F)-(i+p-Kf(/+r) (6.5) 

V+QWV+Q+p) =-(l + Q) 2 (/+®+fi) + (l + Q + K) 2 (/+®) (6.6) 
one gets 



O p 5 

A ong = -— (g',A'| 7AJ |P',A')G A , A ,(t')x 
(2tt) 4 



/ ^ 



, (Q,A|(/+f -ff) 7 "(/+(0)|P,A) 

(/ 2 - m 2 )((; + P-K) 2 - m 2 ){{l + Q) 2 - m e 2 ) + 

(6.7) 



4P S = --^(Q',X'\ 7fl \P',X')G x<>! (t')x 



(Q,A|(/+fy(/+(0+ff)|P,A) 
(I 2 - m 2 )((l + P) 2 - m 2 ){{l + Q + K) 2 ~ m 2 ) + le 

(6.8) 



A l 3 on9 = --^(Q',X'\ 7fl \P',X')G x<>! (t')x 



fi e d n l- 



(Q,X\(/+P-p)rU+®)\P,X) 



(I 2 - m 2 )((l + P-K) 2 - m 2 ){{l + Q) 2 - m 2 ) + ie 

(Q,x\u+?)r(j+m\p^) \ 

J M (l 2 - m 2 )((l + P) 2 - m 2 )((l + Q) 2 - m e 2 ) +iej 



(6.9) 



A l i na = --^(Q',X'\ 7fl \P',X')G x<>! (t')x 



6.2. THE INTERNAL EMISSION GAUGE VARIATION 
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6<n (Q,A|(/+f) 7 "(/+0)|P,A) 



H e d n l 



{I 2 - m 2 ){{l + Pf - m 2 )((l + Qf - m 2 ) + it 
(Q,\\V+P)^V+®+F)\P,\) 



{I 2 - m 2 )((l + P) 2 - m 2 ){{l + Q + K) 2 - m 2 ) + ie 

(6.10) 



Thus it can readily be seen that 

A long + A long + A long + A long = Qj (g n ) 

which proves that {A\ + A2 + A3 + A4} is gauge invariant. The proof for the 
other sums of amplitudes mentioned above proceeds analogously. 

6.2 The Internal Emission Gauge Variation 

It follows from the definition of the vertex function |3.20] and the decomposition 



3.22 that 



ie 5 



167H 

(Q, A| {T^t\m 2 -a,m 2 ) - I*(t, m e 2 , m e 2 )} |P,A) (6.12) 
= ^ T (Q',A'|7M|i" 1 A , )GA^(i')(Q 1 A|7' 1 |i',A>{FF - FF } - 
<Q',A'|^|P',A')G A , A ,(t')(Q,A|^|F,A){FF a + FF b } 



16tt 2 

D 5 



(Q',A'|7jP',A')C? Ai>; (t')(Q,A|7 fl |P,A) x 



16tt 2 

{-4 + 2B (m 2 - a,m ,m e ) - 2B (mo ,m e ,m e ) + 3B (t, m e ,m e ) 

3a 

-3B Q (t',m e ,m e ) - — — {B (m e 2 - a,m Q ,m e ) - B (t' ,m e ,m e )) + 
v + a 

2tC (m 2 1 2m e 2 + t, m 2 , m Q ,m e ,m e ) - 
2t'C${m 2 — a, 2m 2 + t' , m 2 , to , m e , m e ) } — 

" r<Q', A'| ^ |P', \')Gx,x(t r )(Q, A| £ |P, A) x 



16tt 2 

4i' 

{ - -. C (m 2 - a, 2m 2 + t' , m 2 , m , m e , m e ) + 

v + a 

6? - 4a „ . 2 . 10t' , , 

j— - — -^B (m e - a,m ,m e ) - 7- — — r^B {t ,m e ,m e ) 

(v + ay [v + ay 

4 3 

+ — — (- + B (m 2 ,m e ,m e )) } (6.13) 
v + a 2 



The scalar integrals were written in explicit form in 6.13 since the FORM nota- 



tion differs between the form-factor and internal emission routines. It is pivotal 
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for the consistency check that Eq. 6.13 be compared to the FORM output after 



replacing 
Ref. M: 



IfC , which was done in the a3. gauge routine. The result reads, 



167T 



:(Q', A'| 7m \P', \')G\,\>(t')(Q, A| ^ \P, A) x 



It' — a 3t' 
2iCi24 — 2t C134— 2£?i2 + 3B 2 4 + #13 — ti~. — -^34 



t' + a 13 t' + a' 

(Q\ X'\Q \P', X')Gx,x (t')(Q, A| jf |P, A) x 
4 



le 

16tt 2 

4/' 2 

— — 0134-1 

t' + a 

m' 

(f + a)- 



t' + a t' + a 
4q-6f 

B34+ (FT^ Sl3 



12" 



(6.14) 



Eq. 6.14 is expressed with the notation of scalar integrals given in section 
3.3. l| and can be seen to be identical to the result in 6.13] by observing that 
B(j(m^, m e , m e ) = So(m e 2 , too, m e ) — 2. This therefore constitutes the strongest 
confirmation about the internal consistency of the calculation to this point since 
both results were derived in a completely independent manner! 
The relation described by Eq. 5.12] is actually an example of a higher order 
Ward-Takahashi identity jl^, . This can be seen by following Feynman's treat- 
ment of gauge invariance fl29|] , which leads to the following matrix identity 0j for 
the emission of a longitudinal photon from an internal part of the electron line 
with four- momentum P and possible dependence on other external momenta qc. 



K^(P 7 K,q t ) = T(P-K, qi ) - F(P, qi ) (6.15) 
The T-factors for non-zero transfer are directly related to self-energy expressions 



and can be obtained from graphs 3.6, p.4| and 3.12 by removing the virtual 
photon connecting the two fermion lines. Thus, Eqs. 3.12 and 6.15 are two 
equivalent representations of the Ward-Takahashi identity. The Ward-identity 
for zero transfer will now be discussed in the following section. 



6.3 Ward Identity 

Symmetries play important roles in physical theories not only because they 
allow for an elegant and insightful mathematical formulation of physics but also 
because they impose stringent conditions on the possible states of a system or on 
relations between certain quantities of the theory. One important consequence of 
the gauge symmetry in QED is the Ward identity . Since it links self-energy 
and vertex functions in the limit of zero momentum transfer, the calculated 
expressions in the respective chapters should be related according to: 



6.3. WARD IDENTITY 
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1^(0, P 2 ,P 2 ) = - 



<9£(P, mo, / 
dPT 



(6.16) 



The l.h.s. of Eq. |o| can be derived by using the full result of Eq. |3.2l| with 
both fermion legs off-shell and for zero transfer. This was done in the FORM 
routine ffward in the appendix Ref. 124] , and gives: 



r"(0, P 2 ,P 2 ) = {B {P\m ,m e )-B {m \m ei m e )-Am e A C l2 ^-b)r 



(6.17) 



In orde r to obtain the r.h.s. of Eq. 3.16 one has to use the result for X(P, mo,m e ) 
in |3.42l This leads to 



dYi{P, mo, m e 
dP,, 



-(4m e 2 Ci23o + B (m 2 ,m , m e ) - B (P 2 1 m Q , m e ) + 3)7 M 
dB (P 2 ,m ,m e ) 



dP u 



(6.18) 



With 



dB Q (P 2 ,m ,m e 
dP» 

1 d 



in 2 J dPft {I 2 - m 2 )((? + P) 2 - mi) + ie 



Z7T 2 
4 
27T 2 

and using Q 

h = J d n l 
h = [ d n l 



dH 



i» + p» 



(I 2 - m 2 )((Z + P) 2 - m e 2 ) 2 + it 



dH / da: 



x^ + P*) 



o (2 2 + 22Px - m 2 (l - x) + (P 2 - m e 2 )a; + ie) 3 



(6.19) 



1 



(I 2 + 21P + M 2 + ie) x 
{l 2 + 2lP + M 2 + ie) x 



r(A-f) 



i 



r(A) (Af 2 - P 2 + ie) x -% 
P^ I (6.20) 



one finally gets 



(9P (P 2 ,TO ,m e ) 
9P„. 



= 2P^ / dx 



Jo 



P 2 x 2 — (P 2 — m e 2 + m Q 2 )x + m 2 + ie 



48 



CHAPTER 6. GAUGE VARIATIONS 



up to 0{'- 



2P» 



P 2 — m e 2 +m 2 , ?tiq 2 



p2 



Inserting 6.21 into 6.18 and using the relation 
B (rrie,m ,m e ) = 2 + B (mo,m e ,m e ), 



ic 



(6.21) 



(6.22) 



which was already used in section 3.2 



yields the identical result to Eq. 6.17 



This proves that the Ward identity 6.1£ is satisfied by the expressions for the self- 



energy and vertex functions, calculated in this work in completely independent 
ways! 



Chapter 7 

Monte Carlo Results 



7.1 Finite Mass Effects 

In the previous chapters only massless fermions were considered. For high energy 

2 

processes, this is a very good approximation when ^ < 1 as long as the 
radiated photon is not emitted in a direction parallel to one of the fermion 
directions. Within the context of the theory of multiple bremsstrahlung in 
gauge theories at high energies pq] , a method was developed to incorporate 
finite mass effects into a description using massless fermion spinors ||. This 
is especially necessary when running Monte Carlo simulations, where a bulk 
of events is generated in the low angle emission area of phase space and a 
wrong collinear limit would render numerical results meaningless. Following 
the CALKUL collaboration ||, the mass terms important in the collinear limit 
have the following general form in the single bremsstrahlung case: 



da m _ {Q' Q f 



K°dK°dn K dn e 

with 



\A m \ 2 , (7.1) 



KT = -j^f(q-k, Pi ), (7.2) 

where the photon is radiated nearly parallel to q, and / denotes the non radia- 
tive cross section, summed over all polarizations, with the original q replaced by 
q — k. In the case of Bhabha scattering the Born cross section is proportional 
to the following invariant summed matrix element squared j^] : 

/f-+e+ = |>V + « 2 ) (7.3) 

The complete non-radiative cross section for the 0(a 2 ) single bremsstrahlung 
mass corrections is then proportional to 
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f0 _ 
Je~+e+ — 



with 



(l + ^// U )/f- + e + 



(7.4) 



ff°(t) = FF (t)~4irReB(t) 
= -2 + 2Zn^4- 



(7.5) 



From Eq. 7.2 it follows that, when summed over all fermion legs, the finite mass 
terms for the 0(a 2 ) single bremsstrahlung corrections are given by 



\A m i- 



+e+\ 



2e V 
JPKf 



2 r 



4tt 2 



ff°(-2PQ + 2QK) 



[2PP 1 - 2P'K) 2 + (-2PQ 1 + 2Q'K) 2 



{-2PQ + 2QKf 



2e 6 m e < 



e 

47T 



2 ff°(-2PQ-2PK) 



(2QQ' + 2Q'K) 2 + (-2QP' - 2P'K) 



(-2PQ - 2PK) 2 



2e 6 m 2 



(P'K) 2 
(2PP' - 2PK) 2 



1 + -^2ff(- 2P 'Q' + 2 Q' K ) 



[-2QP' + 2QKf 



{-2P'Q> + 2Q'K) 



2e 6 m 2 



An 2 



ff°(-2P'Q' - 2P'K) 



{2QQ' + 2QK) 2 + (-2PQ' - 2PK) 2 
(-2P'Q' - 2P'K) 2 



(7.6) 



7.2 The Soft Limit 

The soft limit of the exact 0(a 2 ) expression is determined by combining all 
helicity- independent (hi) terms of the result 4.3 multiplying the lower order am- 
plitude. Doing this and using Eq. 3.65 and Ref. [G2J gives the following: 



.4 



hi 



16tt 2 

c 2 



- 8m?C% - 2(tCi 24 + t'C&i) + 6B 12 - 6B34} A 

2tt 2 „ -t „ -i! . -t' 



TL 
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Using the virtual YFS-IR function from Eq. 5.11 one gets 



A h Jl - a [Re B(t) + Re B(t')] A 



TL 



I6ir 2 



-4 + Un— % - In— t \ A T , L 
m/ m, 



(7.8) 



In the soft limit t = t' so that 



Af. - 2aRe B(t)A T e L = ^\-l + ln—\A 

2 



hi //"(*) A 



1)1,: 



TL 



\TL 



(7.9) 



where ff° is given by Eq. |7.5|. Taking into account that the cross section is 



proportional to 2Re{A e ~ * A*-}, the soft limit contains the identical form factor 
as the mass correction term in |7.5| . This is another consequence from the YFS- 
theory discussed in chapter |B| which states that in the soft limit, the summed 
and avera ged bremsstrahlung cross section factorizes int o a soft part given by 
S(k), Eq. 5/7, and the lower order cross section, theorem 5.6 . 
It is useful in Eq. 7JS to keep t and t' separately, since t he f ormula holds also 
for large values of k — -g^ 2 — as numerical results in Fig. 7.1 show. In order to 
compare Eq. 7.7 with the result obtained by Fadin et al. ]43(|, one has to use 



the real YFS-IR function given in 5.12 . This then leads to 



A'"- + a 

,,2 



8?r 



B(t, E cut ) + B(t , E cut ) 



e" l —t 

4 + In — =(2lnA - Inx) 



-t' 



m,; 



Al L 



+ln — ^(3 + 2lnA) - AlnA + Inx > A 



\TL 



(7.10) 



with . 



E, 
Ep 



■j^ and A = (in the cm-system Ep = EBeam)- Both Eq. [7^ and 



7.10 are normalized such that they multiply ^da ^ as is the case in Ref. 



The difference to the formula given by Fadin et al. is \ + \ln 2 x. All equations 
are valid for electron-line emission. The other line must be added analogously. 



7.3 Leading Log Comparisons 

At this point all the pieces needed for a Monte Carlo integration are calculated 
and checked on a differential level. It is now necessary to combine the various 
results and include the phase space factors in order to get a normalized differ- 
ential cross section. Following the conventions of Ref. || Q , the expression for 
da o( a+a 2 ) reads . 
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(Q'°) 2 

2 9 lT 5 i 



K°dK°dn K dn e 



7 £ Rz{(A T e L + All )[A e *- + A: + ]} + 



p,\,\' 



(Q'°) 2 



-Y\a 

p,\,X' 



TL 



Alk 



da r 



K n dK°dfl K dn e 



(7.11) 



The amplitudes A e are expressed in Eq. |4.3| ; it is assumed that all hclicity 
degrees of freedom have been summed over and that the initial states were 
averaged. This cross section still contains the IR-mass dive rgen t terms propor- 
tional to In^-, which are identical to those given in Eq. 7/7. Following the 
YFS-prescription of chapter [|, the n orm alized infrared regulated (reg) differen- 
tial cross section for the exact result 4.3 is then given by: 



K°dK°dn K dn e+ K°dK dn K dn ( 



+ 



da° (a) 

^Be{B(t) + B^ } KOdKQdn (7.12) 



An analytic five dimensional integration over 7.12| is extremely involved and 



beyond known calculational techniques. It is also desirable to be able to change 
to a given detector geometry and new experimental cutoff parameters with- 
out having to do a completely new calculation. Both of these difficulties can 
be overcome by using a Monte Carlo (MC) approach [^J. The algorithm im- 
plemented uses the method of weighting Events, i.e. sets of final state 
four-momenta, are generated at random, with a probability for each configu- 



ration to occur given by 7.12. Any experimental situation can then easily be 
simulated by throwing away those events that do not satisfy the experimental 
conditions. The generation of events can be accomplished by first generating 

d& approx 

K u dK"dn K dn e+ ' 



" trial" events according to some approximate cross section K <i fx^dn"^ — ' an< ^ 



assigning to each trial event a weight 



w 



( d 0(a+ a 2 ) \ / d \"1 

I ""reg \ / approx \ „■> 

yK°dK a dn K dsi e + J \K dK°dn K dn e + J ( ' ' 

The exact event distribution is then realized by accepting the trial events with 
a probability proportional to w. From the above considerations it is clear that 
K°dK°dnKdn + mus t be as simple as possible and that the method is only prac- 



tical if the weights |7.13| do not fluctuate too wildly. The approximate cross 
section chosen here is the one suggested by F.A.Berends and R.Kleiss [0 and 
was implemented into Fortran by S.A.Yost. In order to be able to generate the 
phase space variables K° ,9k,^k,0 & + and <i> e + , the integral over da approx has 
to be known and some cuts for these parameters have to be chosen. Here, the 
algorithm of Ref. was used and the cuts correspond to the SICAL W-N 
acceptance H. 
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Figure 7.1: Results of a Monte Carlo integration corresponding to 10 5 events 



reg 



Ofo" 
J reg 



} /&reg 



is a measure for the size of the 



per point; NLL = {a 

subleading and pure (f-) Z terms in the various cases (i). The real photon IR- 
singularity is regulated by a cutoff and cancels out of the ratios. The energy 
cut z m i n is defined in Ref. prl, for example, as is the W-N acceptance. 



The only practical way to calculate J wdo avvrox is to generate uniform distri- 
butions dqi between [0, 1] for each differential factor and write: 



Jdq K odqe K dqg, K dq e ,<ig$ 



(7.14) 



where J is the Jacobian of the transformation 7.14. The integrated cross section 
is then given by summation of all accepted weights and averaging over the total 
number of generated events. 

The results for the integrated exact 0(a 2 ) cross section are presented in Fig. 7.1 
in comparison with several partial results following from the exact expression 
4.3. It can be s een that the ansatz |777| stays within 1%(!) of the exact t-channel 
expression h3 for the chosen cuts. This is a remarkable result given the sim- 
plicity of that equation but should be seen as somewhat accidental, since many 
helicity dependent terms multiplying the tree-level amplitude apparently can- 
cel, without any obvious physical necessity. Taking into account only Tq terms 
is also in excellent agreement with the complete t-channel result and remains 
within 2%. The effect of adding the s-channel, including Z-exchange, can be 
seen to be relat ively small, roughly 5% of the second order single bremsstrahlung 
expression 4.3 for the chosen parameter space. In conclusion, for the accuracy 
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Figure 7.2: A comparison of th e in tegrated cross section with two independent 
calculations, analogous to Fig. 7.1, corresponding to 10 5 events per point. The 



result 4.2 , including s and t channel plus Z-exchange and IR-regularization, is 
seen to deviate in its non leading logarithmic parts from the pure t-channel 
expression given in Ref. |43|| by up to 17%. The LL-ansatz implemented in 
BHLUMI4.01 |39] stays within 8% over the entire range of z min . 



required in the context of 0.1% precision radiative corrections, both t-channel 
approximations 7.7 and using only Tq terms in 4.3 represent perfectly acceptable 
approaches. 

As a further consistency check and also as a measure for the size of the hitherto 
missing subleading terms from the second order single bremsstrahlung calcula- 
tion, the complete result derived in this work with the s-channel and Z-exchange 
is compared in Fig. 7.2 to two independent leading log type results. The size 
of the overall subleading parts is very much dependent of what exactly is in- 
cluded in the leading log ansatz and the authors of Ref. jl3| claim to have all 
terms contributing to a 0.1% precision level in their formula. While an exact 
value for the size of the missing subleading bremsstrahlung parts is difficult to 
specify, Fig. 7.2 suggests that the total contribution of those terms is below 
20% of the total pure second order result. Since the size of the complete second 
order contribution is roughly 0.35% of the Born process for the W-N acceptance 
at SICAL |l5), the subleading terms are found to be small but non negligible. 
The good agreement of the LL-ansatz in BHLUMI4.01, however, suggests that 
the accuracy of the theoretical precision level quoted in Refs. |35[ 
0(0.16%) relative to the Born process, is rather conservative indeed. 
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Chapter 8 

Conclusions 



In this thesis, the sole outstanding contribution to low angle Bhabha scattering 
above the 0.1% precision level, the second order single bremsstrahlung sublead- 
ing terms, have been calculated exactly, i.e. including pure (^) corrections. 
While other leading log and partial results have been published ||, E3I, Eq. 



4.3 represents the first complete, fully differential exact O(a) corrections to 
e~ + e + — > e~ + e + + 7 in the LEP/SLC energy regime. The accuracy of 
the presented result is well below 0.05%, which is the level required in order to 
have meaningful comparisons with measurements of the new luminosity detec- 



tors at LEP |3j . The calculated expression 12 has passed several strong internal 
consistency checks as well as analytical and numerical comparisons with leading 
log type results: 

It was show n in chapter |^ that the expression for the internal emission ampli- 



tudes p.56| is indeed UV-finite and possesses the expected logarithmic virtual 
photon mass dependence as predicted by the Yennie, Frautschi and Suura pro- 
gram. 

Chapter || provided an essential link between technically completely unrelated 
parts of the calculation in deriving the gauge variation for the initial state in- 



ternal emission graph 3.6 by replacing e M — > if M and comparing this with the 
analytical expressions for the on- and off-shell form factors. The demonstrated 
agreement represents a Ward-Takahashi identity for non-zero momentum trans- 
fer. In addition to being a very strong check on the internal emission result, it 
also demonstrates the correct application of the on-shell renormalization scheme 



employed in this work. In section 3.3, the other remaining part of the calcula- 
tion, the self-energy contribution, was linked to the off-shell form factor in the 
limit of zero transfer. The verified Ward identity completes the full circle of all 
the various terms given in this thesis and establishes powerful evidence for the 
overall correctness of the derived expressions. 

In chapter ^ the mass corrections showing a double pole structure, the only 
relevant mass terms in the high energy approximation for a Monte Carlo inte- 
gration, were derived and subsequently added to yield the complete result. The 
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soft limit of Eq. 4.2 was given explicitly in differential form in 7.9 and shown 



to be connected to the lower order cross section according to YFS-theory. By 



combining all hclicity independent terms, an approximate formula 7.7 was found 
that shows agreement with the exact result within 1%(!) after a Monte Carlo 
phase space integration for SICAL-type cuts. It was also shown that the terms 
proportional to the tree-level cross section make up for roughly 98% of the 
second order t-channel contribution and thus contain all relevant leading and 
subleading logarithmic corrections. The size of the s-channel corrections in- 
cluding Z-exchange was found to be as small as expected, around 5% relative 
to the pure t-channel result. Further MC results demonstrated that the exact 
expression 



4.2 relates to other leading log type calculations as follows: 



The IR— regulated ratios {cr l reg — ct^" 2 ) }/<J?e^' 2 \ which are a measure of the 
size of the missing subleading terms in the various partial results (i), revealed 



that the calculation by Fadin et al. lies within 17% of the exact formula 4.3. It 
was furthermore shown that the ansatz used in BHLUMI4.01 |35| is within 8% 
of the answer derived in this work for the specified cuts. Since the size of the 
0(a 2 ) single bremsstrahlung contribution is of the order of 0.35% relative to 
the Born process for cuts chosen to correspond to the experimental luminosity 
detector SICAL at ALEPH, this will have important consequences for the overall 
level of precision of BHLUMI4.X. While several questions about the dependence 
on experimental cuts for the above agreement as well as those relating to the 

it 



technical precision domain remain to be discussed more thoroughly |38|, 37 



is safe to assume that the precision of BHLUMI4.X will be below the 1pm level. 

In light of the results presented in this thesis, it will soon be realistic to achieve 
a theoretical level of accuracy below lpm for the total cross section in low angle 
Bhabha scattering. Together with high precision measurements of second gen- 
eration luminosity detectors at CERN, it will hopefully be possible to determine 
various electroweak parameters in a way that sets stringent conditions on the 
range of validity of this theory and, simultaneously, to open up a new regime in 
the search for physics beyond the Standard Model. 
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